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Studies of ultracold gases in optical lattices^ li 
many disciplines. They allow testing fundamental qus 
turn many-body concepts of condensed-matter physics 
well controllable atomic systems^, e.g., strongly cor: 
lated phases, quantum information processing. Standa 
methods to observe quantum properties of Bose-Einst( 
condensates (BEC) are based on matter-wave interf 
ence between atoms released from trapa 2 -^^^, destrc 
ing the system. Here we propose a new, nondestri 
tive in atom numbers, method based on optical me 
surements, proving that atomic quantum statistics c 
be mapped on transmission spectra of high-Q caviti 
where atoms create a quantum refractive index. This c 
be extremely useful for studying phase transitions 7 -, e 
between Mott insulator and superfluid states, since v; 
ious phases show qualitatively distinct light scatterii 
Joining the paradigms of cavity quantum electrodyna: 
ics (QED) and ultracold gases will enable conceptua 
new investigations of both light and matter at ultinu 
quantum levels. We predict effects accessible in expe 
ments, which only recently became possible 8 . 

All-optical methods to characterize atomic quanti 
statistics were proposed for homogeneous BEG ^' 10111,12 - 



Probe i/ 



and some modified spectral properties induced by BE( 
were attributed to collective emission 9 -^, recoil shifts— 
or local field effects^. 

We show a completely different phenomenon i 
rectly reflecting atom quantum statistics due to sta 
dependent dispersion. More precisely, the dispersion sh 
of a cavity mode depends on the atom number. If t 
atom number in some lattice region fluctuates from : 
alization to realization, the modes get a fluctuating f: 
quency shift. Thus, in the cavity transmission-spectru 
resonances appear at different frequencies directly reflei 
ing the atom number distribution function. Such a me 
surement allows then to calculate atomic statistical qug 
tities, e.g., mean value and variance reflected by spect: 
characteristics such as the central frequency and widt 

Different phases of a degenerate gas possess similar 
mean-field densities but different quantum amplitudes. 
This leads to a superposition of different transmission 
spectra, which e.g. for a superfluid state (SF) con- 
sist of numerous peaks reflecting the discreteness of the 
matter-field. Analogous discrete spectra reversing the 
role of atoms and light, thus reflecting the photon struc- 
ture of electromagnetic fields, were obtained in cavity 
QED with Rydberg atoms^ and solid-state supercon- 
ducting circuits^. A quantum phase transition towards 
a Mott insulator state (MI) is characterized by a reduc- 




FIG. 1: Schematic setup. Atoms are periodically trapped 
in an optical lattice created by laser beams, which are not 
shown in this figure. Additionally, the atoms are illuminated 
by two light modes at the angles #o,i with respect to the 
lattice axis. 



tion of the number of peaks towards a single resonance, 
because atom number fluctuations are significantly sup- 
presse d 17 ' 18 . As our detection scheme is based on nonres- 
onant dispersive interaction independent of a particular 
level structure, it can be also applied to molecule s 19 ! 20 . 

We consider the quantized motion of N two- level atoms 
in a deep periodic optical lattice with M sites formed by 
far off-resonance standing wave laser beamsi. A region 
of K < M sites is coupled to two quantized light modes 
whose geometries (i.e. axis directions or wavelengths) can 
be varied. This is shown in Fig. [T] depicting two cavities 
crossed by a ID string of atoms in equally separated wells 
generated by the lattice lasers (not shown). In practice 
two different modes of the same cavity would do as well. 

As shown in the Methods section, the Heisenberg equa- 
tions for the annihilation operators of two light modes a; 
(I — 0, 1 ) with eigenfrequencies uji and spatial mode func- 
tions ui(r) are 
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ai = -i Ujji + SiD u ^ ai - i5 m D lm a m - na x + T)i(t), (1) 



K 



with Di m = ^ ul(ri)u m (r t )ni, 



2 /A; a , g is the atom-light coupling 
— Lj a are the large cavity-atom de- 



where I ^ m, 5i = 
constant, A; a = loi 

tunings, k is the cavity relaxation rate, r\i (t) — r)ie~ zu ' lpt 
gives the external probe and hi are the atom number 
operators at a site with coordinate rj. We also intro- 
duce the operator of the atom number at illuminated 

sites N K = Yh=i "i- 

In a classical limit, Eq. (fTJ) corresponds to Maxwell's 
equations with the dispersion-induced frequency shifts of 
cavity modes SiDu and the coupling coefficient between 
them SiDiq. For a quantum gas those quantities are 
operators, which will lead to striking results: atom num- 
ber fluctuations will be directly reflected in such mea- 
surable frequency-dependent observables. Thus, cavity 
transmission-spectra will reflect atomic statistics. 

Eq. (TTJ) allows to express the light operators a; as a 
function f(hx, ...,njvf) of atomic occupation number op- 
erators and calculate their expectation values for pre- 
scribed atomic states | *ff) . We start with the well known 
examples of MI and SF states and generalize to any | "J) 
later. 

From the viewpoint of light scattering, the MI state be- 
haves almost classically as, for negligible tunneling, pre- 
cisely (hi)Mi — Qi atoms are well localized at the ith site 
with no number fluctuations. It is represented by a prod- 
uct of Fock states, i.e. |*) MI = FJ^i = \qi, —,Qm), 
with expectation values 



{f(hx,...,h M ))m = fill, —,Qm), 



(2) 



since hi\qi, q M ) = q l \q 1 , ...,q M )- For simplicity we 
consider equal average densities (hi) mi = N/M = n 
((N k )mi = nK = N K ). 

In our second example, SF state, each atom is de- 
localized over all sites leading to local number fluctu- 
ations at a lattice region with K < M sites. Mathe- 
matically it is a superposition of Fock states correspond- 
ing to all possible di stributio ns of N atoms at M sites: 

I*)sf = V.. y/M/M"/^q 1 l...q M l\q 1 ,...,q M ). AI- 

though its average density (hi) sf — N/M is identical 
to a MI, it creates different light transmission spectra. 
Expectation values of light operators can be calculated 
from 



ni, 



■,n M )) S F = ^ w 



Nl 



qil—qia 



f(q 1; ...,3m) ,(3) 



representing a sum of all possible "classical" terms. Thus, 
all these distributions contribute to scattering from a SF, 
which is obviously different from (f(h\, tim))mi © 
with only a single contributing term. 
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FIG. 2: Photon number in a single cavity mode, a, Sin- 
gle Lorentzian for MI (curve A) reflects non-fluctuating atom 
number. Many Lorentzians for SF (curve B) reflect atom 
number fluctuations, which are imprinted on the positions of 
narrow resonances. Here n is smaller than satellite separation 
So (k = 0.15o), N = M = 30, K = 15. b, The same as in a 
but k = So gives smooth broadened contour for SF. c, Spectra 
for SF with N = M = 70 and different number of sites illumi- 
nated K — 10, 35, 68. The transmission spectra have different 
forms, since different atom distribution functions correspond 
to different K. k = 0.05<5n. 



In the simple case of only one mode ag (a\ = 0), 
the stationary solution of Eq. (TTJ for the photon num- 
ber reads 



\Vo\ 



(a p - s a 



(4) 
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where A p = ujq p — ujq is the probe-cavity detuning. We 
present transmission spectra in Fig. [5] for the case, where 
|uo(rj)| 2 = 1, and D 00 = YliLi ^« reduces to Nk- For a 
ID lattice (see Fig. [T]), this occurs for a traveling wave 
at any angle, and standing wave transverse (f9n = tt/2) 
or parallel (9 = 0) to the lattice with atoms trapped at 
field maxima. 

For MI, the averaging of Eq. (TJ| according to Eq. © 

gives the photon number (aJar^Mi, as a function of the 
detuning, as a single Lorentzian described by Eq. (TJ| with 
width k and frequency shift given by <5o(Axi)mi ( ecma l to 
5qNk in Fig. Thus, for MI, the spectrum reproduces 
a simple classical result of a Lorentzian shifted due to 
dispersion. 

In contrast, for a SF, the averaging procedure of Eq. ([3]) 
gives a sum of Lorentzians with different dispersion shifts 
corresponding to all atomic distributions \q±, qx)- So, 
if each Lorentzian is resolved, one can measure a comb- 
like structure by scanning the detuning A p . In Figs.UJa, 
and [2t, different shifts of the Lorentzians correspond to 
different possible atom numbers at K sites (which due 
to atom number fluctuations in SF, can take all values 
0,1,2,. ..,N). The Lorentzians are separated by r5n. Thus, 
we see that atom number fluctuations lead to the fluctu- 
ating mode shift, and hence to multiple resonances in the 
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FIG. 3: Photon number in one of two strongly coupled 
modes, a, Diffraction maximum, doublet for MI (curve A) 
and spectrum with structured right satellite for SF (curve B). 
Structure in the satellite reflects atom number fluctuations 
in SF, while narrow spectrum for MI demonstrates vanishing 
fluctuations. Here ft is smaller than satellite separation 2So 
(ft = O.lrJn), K — 15. b, The same as in a but k = 8q gives 
broadened satellite for SF. c, Diffraction minimum, zero field 
for MI and structured spectrum for SF. Nonzero structured 
spectrum for SF reflects fluctuating difference between atom 
numbers at odd and even sites, which exists even for the whole 
lattice illuminated, K = M. Here k is smaller than satellite 
separation 2r5n (ft = O.lrJn), K = 30. d, The same as in c but 
fc = r5rj gives broadened contour for SF. N = M — 30 in all 
figures. 



spectrum. For larger k the spectrum becomes continuous 
(Fig. [2p), but broader than that for MI. 

Scattering of weak fields does not change the atom 
number distribution. However, as the SF is a superposi- 
tion of different atom numbers in a region with K sites, a 
measurement projects the state into a subspace with fixed 
Nk in this region, and a subsequent measurement on a 
time scale short to tunneling between sites will yield the 
same result. One recovers the full spectrum of Fig. [2] by 
repeating the experiment or with sufficient delay to allow 
for redistribution via tunneling. Such measurements will 
allow a time dependent study of tunneling and buildup 
of long-range order. Alternatively, one can continue mea- 
surements on the reduced subspace after changing a lat- 
tice region or light geometry. 

We now consider two modes with u>q = w\, the probe 
injected only into ao (Fig. [T]) and the mentioned geome- 
tries where D 00 = D n = N K (see Fig. [3]). From Eq. (J), 
the stationary photon number a\a\ = f(n\, tijf) is 



a\a\ 



SfD 10 



Dxo\r} \ 



[A'^-5fD\ D 
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(•5) 



where = A p — #i-Dn. 

In a classical (and MI) case, Eq. ([2]) gives a two-satellite 
contour ([5]) reflecting normal-mode splitting of two oscil- 
lators (ao.i) coupled through atoms. This was recently 
observed 21 for collective strong coupling, i.e., the split- 



ting Si(Dio) exceeding ft. The splitting depends on the 
geometry (see Eq. ([T])) representing diffraction of one 
mode into another. Thus, our results can be treated 
as scattering from a "quantum diffraction grating" gen- 
eralizing Bragg scattering, well-known in different dis- 
ciplines. In diffraction maxima (i.e. u*(rj)itrj(rj) = 1) 
one finds Dig = Nk providing the maximal classi- 
cal splitting. In diffraction minima, one finds D±o = 
J^^ 1 (— l) l+1 ni providing both the classical splitting and 
photon number are almost zero. 

In SF, Eq. ([3J shows that (a{ai)sF is given by a sum 
of all classical terms with all possible normal mode split- 
tings. In a diffraction maximum (Figs. [3^,, b), the right 
satellite is split into components corresponding to all 
possible Nk or extremely broadened. In a minimum 
(Figs. [3fc,d), the splittings are determined by all dif- 
ferences between atom numbers at odd and even sites 



Ef =1 (- 



1) Hi- Note that there is no classical descrip- 



tion of the spectra in a minimum, since here the classical 
field (and {a\a\)yii) are simply zero for any A p . Thus, 
for two cavities coupled at diffraction minimum, the dif- 
ference between the SF and MI states is even more strik- 
ing: one has a structured spectrum instead of zero signal. 
Moreover, the difference between atom numbers at odd 
and even sites fluctuates even for the whole lattice illu- 
minated, giving nontrivial spectra even for K = M. 

In each of the examples in Figs. [2] and [3J the pho- 
ton number depends only on one statistical quantity, 
now called q, f(qi, Qm) = /(<?)• For the single mode 
and two modes in a maximum, q is the atom number 
at K sites. For two modes in a minimum, q is the 
atom number at odd (or even) sites. Therefore, ex- 
pectation values for some state \$>) can be reduced to 
(/)* = S^Lo where py(q) is the distribution 

function of q in this state. 

In high-Q cavities (k <C $o = g 2 /A 0a ), f(q) is given by 
a narrow Lorentzian of width k peaked at some frequency 
proportional to q (q = 0, 1, N). The Lorentzian hight 
is g-independent. Thus, (/) function of A p repre- 

sents a comb of Lorentzians with the amplitudes simply 
proportional to 

This is our central result. It states that the trans- 
mission spectrum of a high-Q cavity (a' a(A p )) directly 
maps the distribution function of ultracold atoms p<sr(q), 
e.g., distribution function of atom number at K sites. 
Various atomic statistical quantities characterizing a par- 
ticular state can be then calculated: mean value (given by 
the spectrum center), variance (determined by the spec- 
tral width) and higher moments. Furthermore, transi- 
tions between different states will be reflected in spectral 
changes. Deviations from idealized MI and SF states^ 
are also measurable. 

For SF, using psf(?) (see Methods), we can write the 
envelopes of the comb of Lorentzians shown in Figs. [2]a,,c 
and[5^,,c. As known, the atom number at K sites fluctu- 
ates in SF with the variance (AA^) 2 = N K {1 - K/M). 
For example, Fig. [^t shows spectra for different lattice 
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regions demonstrating Gaussian and Poissonian distri- 
butions with the spectral width a u — Sq^ (ANk) 2 , di- 
rectly reflecting the atom distribution functions in SF. 
For K w M the spectrum narrows, and, for the whole 
lattice illuminated, shrinks to a single Lorentzian as in 
MI. 

The condition re < do = g 2 /Ao a is already met in 
present experiments. In the recent work&, where setups 
of cavity QED and ultracold gases were joined to probe 
quantum statistics of an atom laser with 87 Rb atoms, 
the parameters are (g, A 0a , k) = 2tt x (10.4, 30, 1.4) MHz. 
The setups of cavity coolin g 23 ! 24 are also very promising. 

For bad cavities (re 3> So = g 2 /Ao a ), the sums can 
be replaced by integrals. The broad spectra in Figs. [2Jd 
and [3p,d are then given by convolutions of p^(q) and 
Lorentzians. For example, curve B in Fig. [2}d represents 
a Voigt contour, well-know in spectroscopy of hot gases. 
Here, the "inhomogeneous broadening" is a striking con- 
tribution of quantum statistics. 

In summary, we exhibited that transmission spectra of 
cavities around a degenerate gas in an optical lattice are 
distinct for different quantum phases of even equal den- 
sities. Similar information is also contained in the field 
amplitudes (00,1) contrasting previous suggestions^ 3 - that 
(ao,i) probes only the average density. This reflects (i) 
the orthogonality of Fock states corresponding to differ- 
ent atom distributions and (ii) the different frequency 
shifts of light fields entangled to those states. In general 
also other optical phenomena and quantities depending 
nonlinearly on atom number operators should similarly 
reflect the underlying quantum statistic s 25 ' 26 ! 27 . 

Methods 

Derivation of Heisenberg equations 

A manybody Hamiltonian for our system presented in 
Fig. [JJ is given by 

H = ^2 fwia}ai+ / d 3 r* t (r) J ff a i*(r), with 

4=0,1 

_ P 2 , , 2 \- <(r)u m (r)a[a m 
Ql ~ 2^T + cl( )+ 9 ^ A ' 

l,m— 0,1 

where ao.i are the annihilation operators of the modes of 
frequencies wo,i, wave vectors k 0j i, and mode functions 
ito,i(r); *(r) is the atom-field operator. In the effective 
single-atom Hamiltonian H a \, p and r are the momentum 
and position operators of an atom of mass m a trapped 
in the classical potential V c \(r), and g is the atom-light 
coupling constant. We consider off-resonant scattering 
where the detunings between fields and atomic transition 
Ai a = — uj a are larger than the spontaneous emission 
rate and Rabi frequencies. Thus, in H a i the adiabatic 
elimination of the upper state, assuming linear dipoles 
with adiabatically following polarization, was used. 



For a one-dimensional lattice with period d and atoms 
trapped at Xj — jd (J = 1,2,...,M) the mode func- 
tions are ito,i(rj) = exp(ijko t i x d + i<f>) for traveling 
and 1*0,1 (rj) = coa(Jko t i x d + (f) standing waves with 
ko.ix — |ko,i j cos#o,i; $0,1 are angles between the mode 
and lattice axes, 4> is some spatial phase shift (cf. Fig.[JJ). 

Assuming the modes ao.i much weaker than the trap- 
ping beam, we expand \&(r) using localized Wannier func- 
tions! corresponding to the potential V c i(r) and keep only 
the lowest vibrational state at each site (we consider 

a quantum degenerate gas): \&(r) = J2i=i biw(r — r i); 
where bi is the annihilation operator of an atom at 
site i at a position r^. Substituting this expansion in 
the Hamiltonian H, one can get a generalized Bosc- 
Hubbard model 7 including light scattering. In contrast to 
"Bragg spectroscopy" , which involves scattering of mat- 
ter waves 4 , and our previous work 2 ^, we neglect lattice 
excitations here and focus on light scattering from atoms 
in some prescribed quantum states. 

Neglecting atomic tunneling, the Hamiltonian reads: 

H=J2 + hg 2 £ (j2 Jftfk) , 

i=0,l i,m=0,l ma Vi=l / 

where hi = b\bi. For deep lattices the coefficients Jjf = 
/ drw 2 {v - Ti)ui (r)u m (r) reduce to J\j = uf (r^Um^i) 
neglecting spreading of atoms, which can be charac- 
terized even by classical scattering 29 . The Heisenberg 
equations obtained from this Hamiltonian are given by 
Eq. (fT]), were we have added a relaxation term. Strictly 
speaking, a Langevin noise term should be also added to 
Eq. (fT]). However, for typical conditions its influence on 
the expectation values of normal ordered field operators 
is negligible (see e.gi^). In this paper, we are interested 
in the number of photons {a\ai) only, which is a normal 
ordered quantity. Thus, one can simply omit the noise 
term in Eq. (TTJ. 

Simple expressions for spectral line shapes in SF 
state 

We will now derive expressions for the spectra pre- 
sented in Figs. [5] and [3] demonstrating relations between 
atomic quantum statistics and the transmission spec- 
tra for the SF state. As has been mentioned in the 
main text, in all examples presented in Figs. [21 and 
the photon number depends only on a single statisti- 
cal quantity, which we denote as q. Using this fact, 
the multinomial distribution in Eq. ([3]) reduces to a 
binomial, which can be directly derived from Eq. ([3]): 
(/>SF = EjL /(«)PSF(g) with p SF (q) - M/[q\(N - 
q)\](Q/M) q {l - Q/M) N ~ q and a single sum instead of 
M ones. Here Q is the number of specified sites: Q is 
equal to K for one mode and two modes in a maximum; 
Q is the number of odd (or even) sites for two modes 
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in a minimum (Q = M/2 for even M). This approach 
can be used for other geometries, e.g., for two modes in 
a minimum and K < Al, where Eq. ([3]) can be reduced 
to a trinomial distribution. 

As a next approximation we consider N, M 3> 1 , but fi- 
nite N/M, leading to the Gaussian distribution psf(<?) = 
l/(y/2na q ) exp [— (q — q) 2 /2a 2 ] with central value q — 
NQ/M and width a q = y/N{Q/M){\ - Q/M). 

In high-Q cavities (k <^ S — g 2 /A 0a ), f(q) is a nar- 
row Lorentzian of width k peaked at some g-dependent 
frequency, now called Ap\ Since the Lorentzian hight is 
g-independent, (f}sF as a function of A p is a comb of 
Lorentzians with the amplitudes proportional to j>sf(<z). 

Using the Gaussian distribution po,-p(q),we can write 
the envelope of such a comb. For a single mode [Fig.[2^,,c, 
Eq. ([!])], we find A? w Soq with the envelope 



(4a (Ap) SF 



aS 



-(A'-A p ) 2 /2^. 



2tt(t u 



where the central frequency A p = SqNk, spectral width 
a u = 8 ^N K {1 - K/M), and a = \r)o\ 2 /n 2 . So, the 
spectrum envelopes in Fig. [2^,c are well described by 
Gaussians of widths strongly depending on K. 

For K — > and K — > M, the binomial distribution 
Psf(q) is well approximated by a Poissonian distribution, 
which is demonstrated in Fig.[2j; for K = 10 and K = 68. 
For K = M the spectrum shrinks to a single Lorenzian, 
since the total atom number at M sites does not fluctu- 
ate. 

In other examples (Figs. [3^ and [3J:), the above ex- 
pression is also valid, although with other parameters. 
For two modes in a diffraction maximum (Fig. [3^), the 
central frequency, separation between Lorentzians and 
width are doubled: A p = 25oNk, A* 2Soq and 

a u = 26 0y /N K (l-K/M); a = \ Vo \ 2 /{2n 2 ). The left 
satellite at A p = has a classical amplitude \r] \ 2 / (4k 2 ) . 
The nonclassical spectrum for two waves in a diffrac- 



tion minimum (Fig. [3];) is centered at A p = SoN, with 
components at A^ « 2SqQ, and is very broad, a u — 

6 VN;a=\ m \ 2 /K 2 . 

For bad cavities (k » So), the sums can be replaced 
by integrals with the same parameters A p and a u as for 
k < So. For a single mode, Fig. EJd represents a Voigt 
contour 



(<4ao(A p ))sF = 



Wo 



2 r°° 



%/27TCr w 7 ( A p - w ) 2 + K 



2 i „2 



For two modes in a diffraction minimum the photon num- 
ber (Fig. Hi) is 



/ t \ W 



where A p = A p — A p , while in a maximum (Fig. [3]b) 

, t v |?7 1 2 Z" 00 uj 2 e-^~ K ^l 2 ^dLo 

(ajai)sF - 



4^2^, Jo [A p (A p -w) + k 2 ] 2 + k 2 w 2 ' 
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